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CURRENT TO A MOVING CYLINDRICAL 
ELECTROSTATIC PROBE 

W. R. Hoegy 
L. E. Wharton* 


ABSTRACT 


The current collection characteristics of a moving cylindrical Langmuir 
probe are evaluated for a range of probe speeds and potentials which are appli- 
cable to earth and planetary measurements. The current expressions derived 
include the cases of the general accelerated current, sheath area limited current 
orbital motion limited current, and retarded current, Eor the orbital motion 
limited current, a simple algebric expression is obtained which includes and 
generalizes the Mott-Smith and Langmuir expressions for both a stationary 
probe and a rapidljr moving probe. For a rapidly moving probe a single formula 
adequately represents both the accelerated and the retarded current. 
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INTRODUCTION 



Cylindrical electrostatic probes have been widely employed on rockets and 
satellites for charged particle temperature and density measurements.^' ^ In 
the future those devices may be used on planetary missions. A common feature 
of the planetary and earth electrostatic probes is the rapid vehicle velocity 
relative to the particle thermal velocity. The ratio of probe speed to thermal 
speed is denoted as the speed ratio. Typical speed ratios for 1000°K ions in the 
mass range 1 to 28 lie between 2 and 10.4 for a probe speed of 8km/eec (earth 
orbit) and between 12 and 65 for a probe speed of 50 km/sec (Jupiter entry probe). 
Speed ratios for 1000°K electrons vary from .046 to .29 for probe speeds between 
8 km/sec and 50 km/sec. 

Mott-Smith and Langmuir^ derived two separate cylinder current expressions, 
one for a stationary probe and the other for a rapidly moving probe. Kanal derived 
some integral formulas valid for arbitrary probe speed and from them, power series 
expressions which are valid only for small speed ratios. Bettinger ^ pointed out 
the need to develop a current expression which fills the void between the formulas 
for the stationary and the rapidly moving probes. In this paper we derive cylinder 
probe current expressions which are valid for a wide range of speed ratios and 
potentials i The current expressions are presented in a form which allow easy 
and rapid numerical evaluation* 



BASIC CURRENT EXPRESSIOI^® 

We review briefly the derivation of the integral expression for the current 
to a moving cylindrical probe* All ear rent formulas are normalized with respect 
to the random probe current and are denoted I, 

I = x/i . (1) 

' random 

where i is the actual current^ and the random probe current is given by; 

X -ANql/iir, (2) 

random ” 2 TT m 

where A is the probe area, N, q, T, and m are respectively the particle density, 
charge, temperature, and mass. 

The probe is assumed to behave as an idealized cylindrical Langmuir probe 
with a coaxial sheath of charged particles which is unperturbed by the relative 
probe-plasma drift motion. It is interesting to point out that the current 
expressions derived here are correct in the limit of high drift velocity - this 
is probably related to the fact that most of the particles collected by the probe 
came from the undistrubed portion of the sheath in the bow region. The probe 
is either guarded or is long enough SO that end effects can be neglected. We 
also assume that no collisions take place within the sheath, and that outside the 
sheath the particles are hot influenced by the presence of the probe. Particles 
which reach the collector surface are assumed to be collected. Under these 
conditions, the cylinder probe current is given by; 
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(3) 
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c c • n f ( c ) , 


where a is the sheath radius, r the probe radius, h a unit normal vector on the 
sheath edge, a the azimuthal angle specifying the location of an area element of 
the cylindrical sheath, and f(c ) is the Maxwell-Boltzmann distribution function in 
a frame fixed on the probe 


f (c) =7r-‘V2 , 


(4) 


where c and s are respectively the particle and drift velocities normalized by 
the most probable particle velocity 


c = V V , s = w V- 

r2kTV r2kT 


m 


( 5 ) 


The dom,ain of integration in Eq. (3) is determined from the condition that the 
particles reach the probe surface using the conservation of energy and angular 
momentum. In terms of the velocity components c ^ , c , c „ , respectively along 
the probe <axis, tangent to the probe axis and the area element at a, and normal 
to the area element at a, the velocity limits are; 


- 00 < c < op 
z 


for retarded 
particles 



(6a) 


/V < c.. < 


- 00 < c . <00 

z 


for accelerated 
particles 




0 < c^< ®, 


where V ^ | e«?^/hT | is the normalized across the sheath potential. The four 
dimensional integral is reduced to a single integral using the change of variable 
c_ = # G 0 s i// , =: ^ s i n and integration over a , , and ijj : 


if d ^ 


(2 g s) 


for retarded particles and: 
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€ d ^ + V I (2 # s) 


/tt Jrv 
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d i X (2 ^ s), 


for accelerated particles. In formulas 7 and 8 s is the component of the speed 


ratio in the plane perpendicular to the probe axis. 
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Two approximations to the accelerating current are oi general interest - 
the sheath area limited current and the orbital motion limited current, The 
current Is approximated by the sheath area limited current when a/r is close 
to unity so that all particles entering the sheath are collected. The current is 
then Independent of the across the sheath potential and depends only on s and a/r. 
To obtain this current the velocity domain (6b) is modified to allow all values 
of the tangential velocity , •” «> < ^ { 
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(9) 


The general accelerated current approaches the orbital motion limited 
current when the sheath is much larger than the probe radius. In this limit, 


a/r ® , the current depends only on V and s; 

^oml =~ f I, (2#s). 

i/n Jo 

The integral current expressions, Eqs. 7-10, are the starting point for the 
derivation of power series representations for the current which are discussed 
in the next section. 

MOVING CYLINDER CURRENT EQUATIONS 

The current collection characteristics for the moving cylinder probe are 
presented in various forms for a range of values of normalized potential V 
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* 


and speed ratio s, The derivations an4 details are found in the appendex 
starting with Eqs. All - A16. 


Sheath Area Limited Current 

The sheath area limited current can be expressed as a single confluent 
hypergeometrlc function j 


Snl 




(11a) 


where 0is defined by Eq, A8. An alternative form has been given by Heatley® 
and *Canal 




a + s') lo 



(11b) 


For small values of the speed ratio, s, the power series repiresentatlon of <!> 
may be used to represent Lai (Eq. A8). For large s, the asymptotic form of 0 
(Eq. A17) gives the representation; 


s large 



( 12 ) 


Hetarded Current 

The retarded current is initially given as a double series in V and s; 
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n^o m»o n ! m j p ^ ^ 


( 13 ) 


Expressing the series in V as a confluent hypergeometrio function and applying 
a Kummer transft.rmatlon (Eq, A21), we obtain the retarded current in the fornx 
of an exponential in V (the retarded current for a stationary probe is e”'^ ) 
multiplied by a power series In the Speed ratio squared: 


'rot 



n - , 


n®o 


v), 


(14a) 


n"0 


(I4b) 


= e”V {l + (l+v)s= + (^-|-l)sn.,.}. 


where (x) is the Laguerre polynomial. Equation 14 is applicable for small 
values of s and arbitrary values of V. This representation of the retarded 
current is particulsyrly useful in describing the retarded electron current to 
satellites and plan^vary probes. 

With the speed ratio, s, kept fixed and small, the retarded current decreases 
exponentially witli increasing potential V| however when the speed ratio is allowed 
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to inotcase, with V kept largo and fixed, the current increases and loses its 
exponential character* The hehavlor is illustrated in the formula of Kanal which 
may bo obtained from Eq. (14a) (see appendix Eqs, A22-A24)i 



As s approaches /V* with V large, the modified Bessel function behaves 
asymptotically as 


e®*^//4 7T s Tv 


so that the exponential behavior is ^ 1 as s - /v. 


Another representation for the retarded current (see Appendix Eq. A25-A29) 
is; 


X 


ret 


\ZWttTT7v 


L 




s /y large, s </v 

vrhere (x) is a polynomial (Eqs. A28 & A29). 

Note that because the potential is retarding, the current is small bat not 
exponentially small as s -♦ </y from below* There is an abrupt change in 
behavior as the boundary s = i^Tls crossed; for s > the drift motion over- 
whelms the retarding potential and the current becomes large as S increases. 
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Again the exponential potential dependence i$ masked by the high drift velocity. 
Tlie expreasion for /V < s is obtained by summing Eq. 13 over s; 



^ 1 * 



(17) 


The fonnulas for evaluating 0 are given in the appendix in Eqs, A30, An 
asymptotic form for s > /Vt s large, is (see Appendix Eqs, A31-A34)j 



s>/V, s large 


(18) 


where (x) is a polynomial of order n (Eqs. A33 and A34). The asymptotic 
behavior of the retarded current is; 

/tt 

Orbital Motion limited Current 


The orbital motion limited currii;:*^ is given by the double series (Eq. All, 
12, 16); 
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The formula suitable for small speed ratios is obtained by summing over 
Index m; 



where ^ is given by Eqs. A7 and A35. The first term in the speed ratio expansion 
(Eq. 20) is the classical orbital motion limited formula of Mott-Smith and 
Langmueir.^ ?or large values of V, the asymptotic series (Eqs. A36-A41) is 
appropriate; 



where the polyonomial Q^, is given by Eqs. A40 and A4X. 

For small values of the potential, another representation is obtained by 
summing Eq. 19 over index rt; 



( 22 ) 
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The first term is identical to the retarded current expression with V replaced 
by -V. The evaluation of - 1/2, 1 ; - s^) is given by Eq, A30 and< 56 (m + 1 , 

1 ; - s*) by Eq. A42. For small values of the speed ratio, s, the retarded current 
and the orbital motion limited current are not identical under the transform 
“ V. However for high s values, the second term in Eq. 22 decays rapidly 
as e“"^ and tlierefore 


s large (V) ^ (- V). 


(23) 


Consequently the asymptotic orbital motion limited eurrent for s ^ > V is obtained 
from Eq. (18) with V - V; 


‘■oml, 



(‘•3 


P /. 1 
n ( 


(24) 


m the appendix it is demonstrated that Eqs. ( 21 ) and (24) are term by term 
identical (Eqs. A43-A45). Therefore these asymptotic expansions may both be 
written as; 


'> (—]' ■ i—f 

wCV + s^)"^ \V + s7 \V+sV 




f 


where 



B(k, {) = (-)' 


4 (- (- 


^ ! - 



2 A 


Writing put the first few terms of Eq. 25 we have; 




/?T 


“S i + 


2 V 


V + s'*' LV + V + S' 


J. 

32 


<V + 


\ / \ + 8 V - 

s^)^ \v + s7 (V + 


8V2 


(V + s2)2j 


(26) 


To first order in v 4 form is identical with the expression; 




/' 




1 S' 


V + V + S- 


(27) 


large or V large 


approximations. For s = 0 we obtain the stationary probe current equation 
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O'"* Eq. 27 
0" 0 


/TTT 


(28) 


and for s > V we obtain the formula; 


lomJkq. 27 


(29) 


r^>v 


which generalizes the Mott-Smith and Langmuir formula for high drift velocity 
by including the random thermal motion in the term 1/2 (the Mott-Smith and 


Langmuir formula is •— /V + , 

fn 


General Accelerated Current 

The power series representation for the accelerated cylinder current is 
obtained from Eqs. A11-A13 but will not be Written explicitly. For small values 
of s, the integral representation is used directly to obtain the representation: 


/tt 


Z] ^ ^ 


). 


(30) 


where f and g^^ are given in the appendix (Eqs. A46-A51). No attempt is made 
to derive an asymptotic form from Eq, 30 for large V since the relation a/r = 
a/r (V) is unknown. For large values of the speed ratio, s, the accelerated 
current reduces to the orbital motion limited current (see appendix Eqs. A52- 
and A53); 
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s large 


(31) 


'acc 


^oml 


Consequently the current to a rapidly moving cylinder probe is given by a single 
formula Eq, 24 or Eq. 25 for both accelerated and retarded particles (with V 
replaced by -V for retarded particles). 


SUMMARY 

The extensive use of the cylindrical electrostatic probe for in situ measure- 
ments of charged particles on rochets ?md earth satjelUtes and its proposed use 
on planetary probes demonstrate the need to understand the effect of the relative 
probe to plasma drift Velocity on the probe volt-ampere characteristics. To 
this end we have derived a number of series representations for the current to 
a moving cylinder probe starting from the four dimensional integral representa- 
tions of the current. Four categories of current have been considered: general 
accelerated current, sheath area limited current, orbital motion limited current, 
and retarded current. Power series in both s, the speed ratio, and V, the 
normalized potential, have been given for all four cases* Simplified representa- 
tions valid for small speed ratio s and arbitrary V well as representations 
valid for large s have been presented, These representations fill the gap between 
the stationary and the rapidly moving probe theories of Mott-Sniith and Langmuir.^ 
One of the more useful of these is the formula for the orbital motion limited 
current: 

- 14 . 


(32) 


oml 


/7 


V 


/ 


V 


V + + 


V + s' 


which is valid for V + > 1. This formula, Eq. 32, can be used to evaluate the 

accelerated ion current for earth and planetary probes. 
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APPENDIX 


Basic Power Series 

We derive a power series i .‘‘presentation for the cylinder prohe current in 
powers of s^ and V from which asymptotic formulas valid for high s or V can 
he oMained, The summation yields series in confluent hypergeometric functions 7 

The evaluation of the cylinder current as given by Eqs. 7-10 veduces to the 
problem of performing the integral: 

X d X /x^ + r Iq (2x s) 

qi/2 

The dependence of I (q, r, s) on s^ is simplified by performing a Laplace trans- 
form in s^: 



This is in the form of a integral representation of the conflttent hypergeometric 

' *4 

fimction^ c; z). Therefore using the relation, 

i/^ (a, c{ z) s z^“® ^ (a - c + 1, 2 ^ c} z), 

we find} 


-2 1 „-pq/i+P A + / i 1. P(q+rA (A5) 

\"T"/ '^\ 2 ' r i + p7 


To facilitate the separation of varis^lps and taking the inverse Laplace 
transform, we partition I into two parts; 

I (q, r, s) = Ij (r, s) + (q, r, s) 

I (q» Tr P) == Ii (f) P) + (q* r. P). 

according to the relation between f and 4>i 


(A6) 


xjj (a, c; z) = <|6 (a, c; z) + z^”® 0 (a - c + 1, 2 - c; z) , 

1 (a - c t 1 ) 1 (a) 


4> (a, c; z) = 


n = 0 


(a)^ z" 
C^)n " ' 


(A8) 


The result of the partition is ; 


18 


M . 



qPf/l+P, 


(A9) 


(q. r, p) - - 


1 (r 4- q)^^^ 
l+p p/5\ 


e"pq/l+P (p 



p^£j^)yA.io) 

I -f p / 


The inverse Laplace transform of I, (r, p) is obtained by expanding it in a 
double series in powers of r and l/p and evaluating the residue at p “ 0. The 
yesult is a power series in s^ and r* 

- r'T (n + m-A) 

I m 1 ^’("^”' 5 ) ^ + 


CD 00 




h»0 m*0 


n 


Similarly I 2 (q> r, p) is expanded in a triple series in powers of (q r), q, and 

l/p. Evaluation of the residue yields: 


q)"» (r + q)"^ P (n + m +1)(A12) 


The four basic currents: accelerated, retarded, sheath area linU^^^ 

orbital motion limited are combinations of I j and 5 
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n ! m t r 


f.i) 


r (m -I- + 1) ^ (n+1) 


lo (q» r, s> = - (r + 


n.m 




I.C. = h <V. h l-^> V. A - 7 Ij / 0, s\ . 


-1 


*• ^1 a2 


- 1 


= (-V. s). 

(A14) 

:..i = 1 1. (0, s). 

(A16) 

I, (V, s) + I, (0, V, s) 

(A16) 


Asymptotic Formulas 


are 


(p (a, c 


2 ) - (a ~ c f l)n T(c) y (a 4 n, z) 

’ ! r (c -a) r(a) ’ 

a - c |: negative integer 

^ (1 - a) r (c) y (c - a + n, z) e‘ 
0 (8j c; z> s 7. - 

^ n I r (a) r (c - a) z<=-“+" 


(A17) 


(A18) 


a ne ga t i v e in t e ge r 


The incomplete gamma funotion, y » may l)e written as 


y (a, X) s a"^ x" («» » + U - x). (A19) 

For sufficiently large z, the following holds, 

r (a) «• 7 (a, z) z^"* e**** (A20) 

The Kummer transform for 4^ is? 


(f> (a, c; z) s e* 0 (c - a, c; - z)« 


(A21) 


Kanal Formula for Ketarded Current 

To obtain the Kanal formula we begin with Eq, (14a), expand and sum the 
series in s: 


Ir.t = 


k"0 


! k r V 2 


cs2\k 

k 


k + 1; - s* , 


(A22) 


Using the Krummer transform we obtain; 


k-O ■ 


Finally we write out the series for 4> and sum over index k: 


/3V 

{jj I 

e-c’*V) 'T 


m"0 


/vj 


(2 s /V) 
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Using the Bessel function representation 




m 


(r. + i, 


2 m + 1 ; 2 X 


). 


(A25) 


We <^tain the formula; 


(I) 




m» 0 


Using the asymptotle expansion of ^ and summing over m, we obtain the result; 


■r«t 


/4 n s /v" 


-(1I-/V)' 


n»* 0 


/ 1 V 

Yfst-v/ {i\ fl^ 

n r “ V24 \2/„ 


The 



2 , i + ni 
2 2 





^ ^2n“3/2 



2 n; 




(A27) 


(A28) 


Where 2^1 * Is a polynomial in s/ /Y” of order minimum of n + 1» 
2m The first three values are: 


% (X) a 1 

(x) S I - 8 X - 8 X* 

<x) 1 - 8 X + 48 X* + 64 x^> 


(A29) 
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Evaluation of (m - 1/2, Ij -s^) 

These functions are defined by the first two along with the recursion relation; 

0 1| . 5“) = e-'Vs (1 + s“) Ij + s2 I, 


0 


(i, 1, - 




0(™ +y li - s“)= _i_ [(2 111 - 52 ) 0^m+i, li-s’j 

m + - 

+ (-2-™) ^(">-5' li -»=)]• 


(A30) 


Asymptotic Form of Ketarded current for s ^ > V 

Substitution of the asymptotic form for 4 > in Eq. 17 yields; 



The sum over index m yields the formula; 


I ^ . — s 

/W 



(A31) 




(A32) 
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A polynomial rfijpresentatlon la obtained \idth the trimsform: 



where P„ (x) - “1/2; x) is a polynomial of order n. The first few 

polynomials are: 



Pj s 1 - 2 X (A34) 

pjsl-8x-8x*. 

Orbital Motion limited Current for Small s 

The confluent hypergeoiiietrio factions t i/' (n - l/2, -1/2; V) are evaluated 
from the first two and the recursion relation: 


"I’ '^) eVerfc(/v) + /v 

, - i { = (1 - 2 v) e'' erfc (A35) 
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■ ■ 


An exact representation for allows us to write the asymptotic form of t/» for 
large Vs 


,h ® "■ [y (a + n, z) 


+ z^“« (-. 2 )" r (c - 1 - n, z) I . 


(A36) 


Consequently we have the asymptotic form: 

/ 1 1 \ R ^)m (“ “ I) 

i/, /n - 1 , - i ; ^ V Y — ^ — — — - 

V 2 2 / Z— t m ! V" 


<A37) 


Substitution into Eq. (20) yields 


R 

=-l Vv ; 

Z— • 



(n + l)n 


("-I 


(A38) 


When s^ > V the series in n converges and may be written as; 




m + 1, m - 1; - -^1 r 


(A39) 


where 2 F 1 has the representation 


,F,= 1 + 


1/ 2- 2m 


3 , s2\ 


where Q„(x) = 2 Ft •" hi; 1; x) is a polynomial of order m. The first 

three terms are given as; 



Q2*= 1 + X - 


Evaluation of 4>{m + 1, 1; -s*) 

Using Eq. A21 we find: 

(m + 1 , 1; - s^) - e”“ cf) (- m; 1; s^) 
= e“»' L<») (s^V, 

where Lj,°^x) is the Laguferre polynomial. 


(A41) 


(A42) 


Term by Term Equivalence of Expansions 

The term by term equivalence of the two orbital motion limited series for 
V > s ^ and s ^ > V follow from the proof of the equality; 


The proof is accomplished by rearranging the polynoniiais. 
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(-r Q. 


k “ 0 


» (-"’)k(|-"') (- . 

= (s 2 )" ' S' 

k ! k I 


k ! k ! 


Using the relation; 


(a - m)„ s ~ (-) . 


. » (D™ f--^) (- ij 


ksQ k 1 k I (1) 


.-K (- s) 


Changing the dummy index m - k-* k, and using 


(m - k) I 


= (-)" (- mX 


w© find 


= (s^)" 


(~l 

^ kiY-i) 

/ \ V 2 ^ 

-T 

V •6/m _ / V 4 


= P /- 

^ ■ m t "> I 
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Accelerated Current ITunctions 


The functions and used In the '>jpeed ratio expansion for are given 
by; 


(1 + 7^) V erfc / (1 + 


fl =^Y+V r^) /(l +7^) V 


+ - V j e*Y er fc /(I + ^> V, 




r 


2 - 


- 1 


fnt2= ('’+!-■'') f„+l +V(n + 1) f„ 


+ y'^ ((1 + 7^) Y^n 


(A49) 


grt = 


/F 


erf (7 /V) - 7 /V* 


(A50) 




^ ^^2 Y^ n+3/ 2 g -r V 


Form of for Large s 

From Eqs. A13 and A16 we note that 

= ^oml 4 algebric sum of three 12 's 

The asymptotic behavior of 1 2 (q, r, s) for s large and fixed q and r is obtained 
by summing Eq^ A12 over index n: 

I 2 (q, r, s) = - (r + q)’''2 ^ ^.5}" S'' . ?6 (m + -t + 1, 1; _ s’), 

m ! r W + ij 

and 

0 (m + + 1, 1; - s^) = e“® <p (- m - 't * 1; + s^) 

which decays exponentially for large s. 

The hypergeometric functions in Eqs. A27, A32, and A39 were expressed 
as polynomials using the theorem: 

F(a, b; c; x) = (1 - F(c - a, c - b; c; x). 

An alternative method of evaluating these functions is directly through their 
recursion relations. The recursion relations may be obtained either from the 
contiguous function relations or by the use of successive raising and/or lowering 
operators. In the latter method, any two of the hypergeometric functions are 
expressed as products of raising and lowering operators acting on the third 


function* Derivatives of order two’ or more are written as derivatives of order 
1 and 0 using the differential equation* Finally tlie recursion relation is obtained 
by eliminating the first derivative. iiUthe^ two resulting equations. 


